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ABSTRACT 

It is estimated that the Weibel instability is not generally an effective mech¬ 
anism for generating ultrarelativistic astrophysical shocks. Even if the upstream 
magnetic held is as low as in the interstellar medium, the shock is mediated not 
by the Weibel instability but by the Larmor rotation of protons in the back¬ 
ground magnetic field. Future simulations should be able to verify or falsify our 
conclusion. 

Subject headings: instabilities-magnetic fields -plasmas -shock waves-gamma 
rays: bursts 


1. Introduction 

There is large literature on gamma-ray burst (GRB) afterglows based on the assumption 
that the X-ray, optical and radio afterglows are the synchrotron emission from relativistic 
electrons Fermi accelerated at the forward shock of the blast wave (see recent reviews by 
Meszaros 2002; Zhang Sz Meszaros 2004). A longstanding difficulty with this assumption 
has been the inferred magnetic field needed to fit the afterglow data typically requires that 
the magnetic energy density exceeds by many orders of magnitude that which would be 
expected from the shock compression of the interstellar magnetic field of the host galaxy 
(Gruzinov & Waxman 1999; Gruzinov 2001). Many authors have therefore assumed that 
the shock somehow manufactures field energy to meet this requirement, but no convincing 
mechanism has been proposed to date. One mechanism discussed is the Weibel instability 
(Medvedev & Loeb 1999, Silva et al. 2003; Frederiksen et ah 2004; Jaroschek et ah 2005; 
Medvedev et al. 2005; Kato 2005; Nishikawa et al. 2005), which has the fastest growth rate 
and produces relatively strong small scale magnetic field even in an initially non-magnetized 
plasma. It is expected that the thermalization of the upstream flow could occur via scattering 
of particles on the magnetic fluctuations. In the electron-positron plasma, the instability 
does generate the magnetic field at about 10% of the equipartition level and does provide 
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the shock transition at the scale of a dozen of electron inertial length (Spitkovsky 2005). 
However, simnlations of colliding electron-proton flows show that while the electrons are 
readily isotropized, the protons acqnire only small scattering in angles after passing the 
simnlation box (Frederiksen et ah 2004). How long the held persists after the shock is also 
an important qnestion (Gruzinov 2001) bnt here we discuss whether the Weibel instability 
can even cause the shock in the hrst place. 

Moiseev & Sagdeev (1963, see also Sagdeev 1966) analyzed the structure of the nonrela- 
tivistic Weibel driven shock and found that the width of the shock transition should be very 
large because electrons easily screen proton currents thus suppressing development of the 
instability. Failure of the Weibel instability to preempt other shock mechanisms, except for 
very large Alfven Mach numbers, has been discussed in the context of nonrelativistic shocks 
by Blandford and Eichler (1987). Here we estimate the width of the Weibel driven shock in 
the ultra-relativistic electron-proton plasma. Although based on a number of physical as¬ 
sumptions about the behavior of the plasma parameters at the non-linear stage of the Weibel 
instability, such analytical scalings are necessary in any case, because, by evident reasons, 
simulations of plasmas are possible only with artihcially low proton-to-electron mass ratios 
(e.g., Frederiksen et al. (2004) took nip/me = 16). In this paper, we present the parameters 
in physically motivated dimensionless form and we believe that our assumptions could be 
checked by numerical simulations. Only by combining numerical simulations with analytical 
scalings can we achieve reliable conclusions about the properties of real shocks. 

As a model for the shock formation, we consider collision of two oppositely directed 
plasma flows. Eventually two diverging shocks should be formed with plasma at rest between 
them. However at the initial stage, the two flows interpenetrate each other exciting turbulent 
electro-magnetic fields. Particles are eventually thermalized by scattering off these turbulent 
helds. As electrons are thermalized relatively rapidly, we consider development of the Weibel 
instability in two proton beams propagating through relativistically hot isotropic electron 
gas. We estimate the proton isotropization length in such a system and conclude that the 
Weibel-mediated shocks are so wide that even in the interstellar medium, the shock should 
be formed at the scale of the Larmor radius of the proton in the background magnetic held. 

The article is organized as follows. In sect.2, we hnd the growth rate of the proton 
Weibel instability. Saturation of the instability is considered in sect.3. In sect.4, we exploit 
the obtained results in order to estimate the width of the Weibel-mediated shock transition. 
Section 5 contains the discussion. 
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2. The Weibel instability 


Let us consider two proton beams of equal strength propagating in opposite directions 
along the 2 :-axis. For the sake of simplicity, let us adopt the waterbag distribution function: 

1 


f;(p) = 


2tip 


2 [Kpz - p\\o) + sipz + piio)]0(p1o - pi); 
±0 


( 1 ) 


where Q{x) is the Heaviside step function, p± = a/Px + Py- Assume that the beams propagate 
through an isotropic relativistically hot electron plasma with the distribution function Fe{p). 

This conhguration is known to be unstable because a small transverse magnetic held 
B = Hxexp {—tu + iky) would drive the oppositely moving protons into current layers of 
opposite sign, which reinforce the initial held (see, e.g., Fig.l in Medvedev & Loeb 1999). 
Evolution of the electromagnetic held is governed by Maxwell’s equations 


_ _ <9B _ „ (9E ^ . 

VxE = -^i VxB = ^ + teji 

which are written in Fourier space as 

kE = uB] i{ujE — kB) = Anj. 
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Note that only components of E and j are present in this conhguration. The current 
density, j, may be found from a solution to the linearized Vlasov equation 


i{(jj — kVy)SEp^e = ±e 
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As usual (e.g., Krall & Trivelpiece 1973), the condition for the set of equations (3-5) to have 
a nonzero solution yields the dispersion equation 
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is the susceptibility of the plasma species a. 

Substituting the proton distribution function (1) into Eq.(7) yields 
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where ilpp = Aire^/mp'j is the relativistic proton plasma frequency, 7 = + {p]_Q + p|g)/mp. 

Below only the strongly anisotropic, highly relativistic case is considered, P||, 7 3> 1. 

Then one can neglect the hrst term in the curly brackets. 



The susceptibility of isotropic electrons is written as 
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Here we used the spherical coordinates in the momentum space and the Plemeli formula. The 
Weibel instability operates in the low-frequency limit, oj <^k. In this case the imaginary part 
of Xe dominates because at cu —0, the principal value of the integral in x goes to zero. The 
physical reason is that only electrons with small Vy contribute to the susceptibility because 
other electrons ’’see” rapidly oscillating held as they move in y direction over a distance 
larger than 1/k for the time ~ 1/oj. Now one can write 


( 10 ) 

where Qy = Ane'^n/T, 1/T = Stt f F^pdp. The parameter T is equal to the electron tem¬ 
perature if the electron spectrum is Maxwellian and T 3> mg. Note that in their analysis of 
the proton Webel instability, Wiersma & Achterberg (2004) erroneously used the expression 
Xe = ~^pel^‘^i which is valid only in the high frequency limit, ct; 3> /c, and therefore is 
irrelevant to the case of interest. 


Now one can write the dispersion equation ( 6 ) in the low-frequency limit: 
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In the limit of negligible angular spread of the proton beams, p±o 
equation is reduced to a simple cubic equation 


( 11 ) 

0 , the dispersion 
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where a = 7 rr 2 pp/ 4 r 2 pg = 7 rmp 7 / 4 T, k = k/Qpp, x = —zu/flpp. The system is unstable 
provided > 0. In the small- and long-wavelength limits, the growth rate 77 = = xflpp 


is found as 


f ka k 

\ Vtpp] k > a^k^Vtpp. 


(13) 


The full solution to Eq.(12) is shown in Fig.l. One can see that if T ~ yuie, as one 
can expect within the shock structure, the most unstable are short-wave perturbations, 
k > {mp/rrieY/^klpp. 


When the angular spread of the beams increases, the growth rate of the instability 
decreases (Fig.2). The threshold of the instability may be easily found by substituting cu = 0 
into Fq.(ll); this yields 

P-Lo _ V^klpp 
P||o k 

So at small wavelengths, where the growth rate is maximal, the instability stops after a small 
spread in the angular velocity distribution is achieved. 


3. Stabilization of the Weibel instability 


The dispersion relation ( 6 ) was found assuming that the particles are nonmagnetized, 
i.e. their trajectories are nearly straight. The instability saturates when this condition is 
violated either for protons or for electrons. As a result of the instability, current hlaments 
are formed along the direction of the proton motion. The magnetic field forms a sort of 
cocoon around these hlaments and eventually traps protons within the hlaments; then the 
instability stops (Yang, Arons & Langdon 1994; Wiersma & Achterberg 2004). The quiver 
motion of a proton within the current hlament may be described by the linearized equation 

= eu25(^,t); (15) 

where ^ is the proton displacement in the transverse direction. Near the axis of the hlament, 
the magnetic held may be written as i? ~ k^Bo, where Bq oc exp? 7 t is the amplitude of 
the perturbation. Then Fq.(15) describes oscillations in the transverse direction with the 
frequency Uq = y/eBok/rrip'y and the (growing) amplitude ^ = cBq/{'^ mprY)- The proton is 
trapped within the hlament when the oscillation amplitude becomes less than 1/k or, which 
is the same, frequency cuq exceeds the growth rate of the instability, rj. This occurs when the 
magnetic held reaches the value 
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One can consider electrons as non-magnetized if their Larmor radius exceeds the char¬ 
acteristic scale of the unstable perturbation (Moiseev & Sagdeev 1962; Sagdeev 1966). This 
condition is violated when the held reaches the value 


a. ^ 


( 17 ) 


Then the magnetic held becomes frozen into the electrons and the magnetic hux does not 
grow more. 

The obtained limits on the magnetic held are shown in Fig.3. The instability develops 
until the magnetic held Bgat = niin(i?fr, i^trap) is reached. One can see from Fig.3 that the 
maximal held is achieved in perturbations with the wave number 


ko — O' ^ ^pp — \l A ^pe 7 


(18) 


the corresponding wavelength is of the order of the inertial length of electrons. The energy 
of the generated held scales as the energy of electrons: 




1 / mp'jQ. 


p'“pp 1 = ZpiT. 

n 


(19) 
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Within the shock structure, one can conveniently normalize the electron temperature as 

T = rme7/3; (20) 


where r is a dimensionless parameter; r = 1 means that the average electron energy remains 
the same as upstream of the shock. Now one can estimate the fraction of the upstream 
kinetic energy transformed into the energy of the magnetic held as 

2rme 

€b — -• 

STinip 



4. The width of the shock wave 

When two oppositely directed plasma streams collide, the Weibel instability generates 
small scale magnetic helds; particle scattering oh these magnetic huctuations provides an 
isotropization mechanism necessary for the shock transition to form. The electron streaming 
is halted easily whereas protons still plow on through an isotropic electron gas. The shock 
is formed on the scale dehned by slow dihusion of protons in the momentum space (Moiseev 
& Sagdeev 1962; Sagdeev 1966). 
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Fig. 3.— Saturation magnetic field as a function of the wave number. The limit i?trap due 
to the proton trapping is shown by solid line; the limit Bit- when the field becomes frozen 
into electrons is shown by dashed line 
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Let us denote the transverse scale of the magnetic inhomogeneities by d and the ampli¬ 
tude of magnetic fluctuations by B. The proton is scattered over a characteristic correlation 
length by the angle 


66 


eBd 

sin OrUp'-f ’ 


( 22 ) 


where 6 is the angle the proton makes with the flow direction. Here we take into account 
that the Weibel instability generates strongly elongated current hlaments so that the proton 
passes the distance / = d/ sin0 within the same hlament (of course I and 66 remain hnite at 
6^0 however we are interested in isotropization scale, which is determined by 6* ~ 1). The 
angular diffusion coefficient is estimated as 


^ _ e^B^d 

I sin dnip'y'^ ’ 

which yields the isotropization scale 

_ 1 

d\eB ) ' 


(23) 


(24) 


Motivated by the estimates expressed in Eqs. (18), (19), and (20), we normalize the charac¬ 
teristic inhomogeneity scale, d, and turbulent magnetic held amplitude, H, as 


B = rnic'i/3-, d = 


(25) 


where r, and C, are dimensionless parameters. Now the shock width may be expressed as 


/ STTrUp \ 1 

\4Tme) 


(26) 


The estimated shock width, L, is thus very large compared to the proton inertial length 
1/Hpp assuming that the electron temperature as well as the scale and amplitude of the 
generated magnetic held do not exceed signihcantly their hducial values, i.e. r ~ ~ C ~ 1- 

We now explain why this is expected to be the case. 

If r ~ mp/rrie 3> 1, then, although the shock width could then be brought down to 
the proton skin depth, this would beg the question of how the electrons are heated, which 
is merely passing along the question of a shock mechanism. Similarly, we expect ^ 1 
because we know of no reason to expect that magnetic held would grow above the limit (19). 
On the contrary, such a small-scale magnetic held should decay (Gruzinov 2001). There 
are evidences for hierarchical merging of current hlaments (Silva et al. 2003; Frederiksen 
et al. 2004; Medvedev et al 2005; Kato 2005) so that ( might exceed unity. However, the 
merging should be accompanied by the held dissipation because the magnetic held between 
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two parallel adjacent currents change sign and therefore should dissipate when the currents 
merge. So possible increase of ( would be compensated by decrease of In any case, at 
the scale larger than the electron Larmor radius, which is of the order of l/(Q;^/^r2pp), the 
held is already frozen into the electron gas therefore ( could hardly ever grow signihcantly 
beyond unity. One should also take into account that the current hlaments are unstable to a 
kink-like mode (Milosavljevic & Nakar 2005), which also stimulates the held decay. Therefore 
we believe that there is no reason for the Weibel-driven shock transition to be signihcantly 
narrower than Eq.(26) predicts. There is, however, reason to suspect that the transition is 
even more gradual than predicted by Eq.(26); this is decay of the small scale magnetic held. 

The highest resolution published simulations of shell collisions with mp/me = 16 (Fred- 
eriksen et al. 2004) do show that while electrons are readily isotropized, the proton beams 
achieve only a small angular spread when passing the simulation box of the length dO/hlpp. 
These simulations also show that the average electron energy does not grow so that our 
conjecture r ~ 1 seems to be correct. The experiment duration, 120/12^^, was three times 
larger than the particle crossing time; by the end of the simulations the spatial wavelength of 
the magnetic fluctuations achieved one half of the width of the simulation box. In physical 
units, this is written as A = 20/Qpe-, k = O.Slflpe so that C, = 3. Although one can not 
exclude that the pattern growth was frustrated by the hnite size of the simulation box, we 
believe, by the reasons outlined above, that C, would not grow signihcantly in any case. 

According to the estimate (26), the full shock transition is too wide to be simulated 
numerically even with a moderate mass ratio mp/me > 10. On the other hand, the scaling 
(26) may hardly be ever applied to the case mp/me ^ 10 because it was obtained under the 
assumption that the proton and electron scales are well separated, i.e. that mp/me ^ 1- 
Therefore direct numerical check of this scaling is very difficult. Nevertheless it would be 
very useful to follow behavior of the parameters r, ^ and ( in numerical simulations even 
with a low mass ratio. Even 2.5D simulations of the proton Weibel instability in the isotropic 
electron gas would clarify the behavior of these parameters in highly nonlinear regime. 


5. Discussion 

The estimate (26) was obtained under the assumption that there is no magnetic held in 
the upstream how. If the how is magnetized, a shock transition may be formed at the scale 
of the proton Larmor radius; therefore the above estimates are valid only if eBoL/mp'-f < 1, 
where Bq is the magnetic held in the upstream how. One can conveniently characterize the 
magnetization of the how by the parameter a = Bq / {A'nmp'yn). Making use of Eq.(26), one 



-11 - 


finds that the shock may be driven by the Weibel instability if 

= 1.5x10-“5VC=, (27) 

\ STirrip J 

If the shock propagates through the interstellar medium, the magnetization exceeds the right- 
hand side of Eq.(27) by factor about 30: a = 5 ■ 5 n“y where B = 10“^'®-B_5.5 G is 

the interstellar magnetic held, n cm“^ the gas number density. This suggests that forward 
shocks that presumably produce GRB afterglows are mediated not by the Weibel instability 
but by the Larmor rotation of protons in the background magnetic held. This does not mean 
that the Weibel instability does not work at all. On the contrary, it does develop and may 
well create some small scale magnetic held that is much stronger than the background held. 
The scattering oh these magnetic huctuations results in dihusion of the protons in angles. 
However, because the scattering does not manage to isotropize the protons at the scale less 
than the Larmor radius it is hard to see how the Weibel instability could convert the kinetic 
energy to another forms. The strong dependence of the estimate (21) on the parameters 
r, and C, makes accurate determination of their values, presumably by simulations, crucial 
to solidify this conclusion. 

The estimate (27) shows that a fraction ~ 10“^ of the total energy is converted into 
the magnetic energy unless the electrons are heated additionally within the shock structure. 
The generated small-scale held should decay (Gruzinov 2001; Milosavlevic & Nakar 2005) so 
that es may be even lower. According to Panaitescu & Kumar (2002) and Yost et ah (2003), 
the observed spectra and light curves of the GRB afterglows imply ~ 10“^ -i-10“^ in most 
cases. Eichler & Waxman (2005) demonstrated that the above estimates may be rescaled 
such that the observations are htted with values of es that are smaller by an arbitrary factor 
/, nie/mp < / < 1. Taking this into account one can see that the Weibel instability could 
provide the necessary held unless the held decay is too strong. On the other hand, as the 
global structure of the shock transition is dictated by the Larmor rotation of the protons in 
the background held, some new physics could come into play. 

A presumably important physical mechanism is the synchrotron maser instability at the 
shock front. This instability generates semi-coherent, low-frequency electromagnetic waves 
(Gallant et al. 1992; Hoshino et ah 1992; Lyubarsky in preparation). In low magnetized 
hows, the amplitude of these waves exceeds the amplitude of the shock compressed back¬ 
ground held. In this case, relativistic particles radiate in the held of the waves via nonlinear 
Gompton scattering (e.g., Melrose 1980, pp. 136-141). The power and characteristic frequen¬ 
cies of this emission are similar to those for synchrotron emission in the magnetic held of the 
same strength therefore the observed gamma ray bursts afterglows may be attributed to the 
nonlinear Gompton scattering oh the electromagnetic waves generated by the synchrotron 
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maser instability at the shock front. It is beyond the scope of this paper to redo afterglow 
theory with the spectrnm of low freqnency electromagnetic waves that is expected from this 
instability. 
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